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Abstract 

We have previously shown an analysis of our dimer model in the over-damped regime to show 
directed transport in equilibrium. Here we analyze the full model with inertial terms present to 
establish the same result. First we derive the Fokker-Planck equation for the system following a 
Galilean transformation to show that a uniformly translating equilibrium distribution is possible. 
Then, we find out the velocity selection for the centre of mass motion using that distribution on 
our model. We suggest generalization of our calculations for soft collision potentials and indicate 
to interesting situation with possibility of oscillatory non-equilibrium state within equilibrium. 
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In a previous communication [1] we have shown a simplified calculation at the over- 
damped limit of a dimer model with particles undergoing hard collisions. With this model, 
we demonstrated the possibility of achieving sustained broken symmetry for the system of 
dimer in contact with a heat bath in the absence of any external drive. As a result of that 
sustained broken symmetry under equilibrium conditions, the dimer system shows directed 
transport. We have already given a lengthy discussion in ref.[l| on why equilibrium does 
not forbid possibility of directed uniform motion, how the model be generalized for soft 
collisions and the possible role of dissipation at collision, when the dimer-bath system as a 
whole remains isolated from surroundings, causing uniform directed motion in equilibrium. 
In this paper, we do a detailed calculation treating the full model (containing inertial terms) 
to show that our previous results at over-damped limit was correct. Our present results, 
being identical to the over-damped ones we have already shown, will establish that the 
hardcore collisions we had considered in the previous case [l| through boundary conditions 
was a correct way of dealing with it at the over-damped limit. 

Our dimer model [1] with hardcore collisions between particles of the dimer would look 
like the following. 
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Z — V\ — v 2 



In the above model z = x\ — x 2 where x\ and x 2 are the positions of the particles of the 
dimer such that X\ > x 2 . As we have shown in the previous article [l| that, the requirement 
of having two different damping during the particles of the dimer are on flight (i.e. when the 
harmonic force is working on them) and during the period they physically collide with each 
other is the basis of the sustained symmetry breaking. Here, since we are going to analyze 
the system with hardcore collisions, for the sake of simplicity, we consider that the time over 
which actual collision takes place is negligibly small. For example, the time of collision is 
of the order of time of a single collision of the dimer particles with a bath molecule (say). 
Within this consideration, we can safely drop the noise and the damping during collision and 



we also drop the harmonic interaction during collision because this force is proportional to 
the separation of particles which is practically zero during collision, we emphasize on the fact 
that, taking this hardcore limit is for the sake of simplicity and one can readily generalize 
this model for soft-core collisions where one must take a different damping constant during 
the time the collision happens and the strength of the noise in accordance with this damping 
during collision to have the symmetry breaking. Having the broken symmetry there, one 
would find the qualitative result (i.e. directed transport in equilibrium) is present, of course, 
with some more variety. We would talk about this generalization in the end of this article. 
Note that, < r)i(t) >= and < T)i(t)r)j(t) >= 5ij8(ti—t 2 ). The mass of the particles are taken 
the same and unity for the sake of simplicity which would make the numerical simulation 
of the above model implementing hardcore collisions straightforward. In the above model, 
consider that the parameter a in the collision potential is so large that the particles have no 
chance to pass through each other for a given k B T . 

Following a Galilean transformation to a moving frame moving with a velocity V cm one 
can rewrite the above model in the form 



lim9(Z - e)[ - (1 - P){v x - V cm ) + ^2(1 - f3)h B T Vl (t)} 
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v 2 = \im&(Z - e)[ - (v 2 - V cm ) + ^2k B Tr] 2 {t)\ 

e— >0 

e(z-e)^ + (l-6(z-e))(a-^)] (4) 
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Considering the time dependent probability of a state (z,v\ } v 2 ) being given as 

P(z, Vl , v 2 , t) =< 8(z - z(t))5( Vl - v 1 {t))8{v 2 - v 2 (t)) >, (6) 
where <> means average over noise. In general, for the Gaussian noise rj, the probability 
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distribution of the noise is given by P[rf\ = Cexp(— J dt ^-f? T ), where C is normalization 

to 

constant and T being the damping is to be taken (1 — f3) for rji and unity for rj 2 in our actual 
calculations. Its generally taken that the time gap tf—t = r ^> r c , where r c is the correlation 
time for noise which is practically zero for the (^-correlated noise approximation. The noise 
average of a quantity X[rf\ is obtained by evaluating the functional integral f T>[t]]X[t]]P[t]], 
where the discrete time steps are r. Since, in our hardcore collision limit, we are considering 



the collision time is of the order of r c , it would not practically affect the calculations and a 
very standard result would follow as one does the relevant averaging over the noise. Taking 
into consideration the continuity equation for the normalized probability P(z,v 1 ,v 2 ,t), we 
can easily derive the Fokker-Planck (FP) equation for the above mentioned model on a 
co-moving frame from 

dP d 

— = -— < S(z - z(t))<S(ui - ui(*))fl(u 2 - v 2 (t))(v 1 - v 2 ) > 
° < S(z - z(t))S( Vl - ui(t))<J(u 2 - v 2 (t)) 
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The FP equation that follows would now look like 

dP , ,dP d$(z) (OP dP\ 

= lime(Z -e)(l-0)± ( {V1 - Vcm )P + k B T d £ 

+ lime(Z - e) A f ( U2 - \/ cm )P + ^t|^) . (8) 

The form of the above mentioned FP equation is independent of the choice of the internal 
field potential ®(z). So, long as the field is internal, the FP equation can be written in 
the above mentioned form because one can always write down the model in terms of the 
internal coordinates and the individual particle velocities. Since the system would settle 
in equilibrium at the minimum of the potential (of course with fluctuations present), &(z) 
being internal the system and the field can move together and the above mentioned form of 
the FP equation simply mentions this fact without giving a selection for the velocity of CM 
which has to be found out from the dynamics. 

The steady state or equilibrium solution of the above FP equation is the standard one 
but with a constant drift and can be written down as 

P(Z,v 1 ,v 2 ) = Ne~^ (z)/kBT e~ {vi ~ Vcm)2/2hBT e~ {v2 ~ Vcm)2/2kBT (9) 
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where N is the normalization constant and $(z) = ]im € -^ [Q(z — €) 9 ^- + (l — Q(z — e))(a — 9 f)]. 

having found the equilibrium solution of the relevant FP equation to our model, we 
would now use this distribution to find the various average quantities to get a selection of 
the uniform average drift V cm . To do that, we would first write down our model into the 
form involving V cm = {v\ + i>2)/2 and V z — V\ — v 2 . So, now the model looks like 

V cm = lim Q{Z - e)[(| - l)V cm + jV z + U] (10) 

8 az 2 
V z = lim Q(Z - e) [(£ - 1)V Z + (3V cm + £ g - 2az] - 2 lim — 5(z - e) 

-2 lim [-(1 - G(z - e)) + (a - — - e))l (11) 
i = K (12) 

In the above form of the model £ cm = (^2(1 - pJk^Trj^t) + ^/2k^Tr] 2 {t)) /2 and ^ = 



a/2(1 - /3)k B Trji(t) - v / 2k B Tr] 2 (t). Its easy to see that < £ cm >=< ^ 2 >= 0, but, £ cm and 
^ 2 would have cross correlations and that is why we have taken the other form of the model 
involving vi and v 2 to derive the FP equation. Now the analysis is straight forward, replace 
for Q(z — e)V z in Eq.ll from Eq.10 then take the average with respect to the above mentioned 
probability distribution. The < V cm > now comes out the constant V cm which is the uniform 
drift of the frame. Set all the average accelerations equal to zero (equilibrium) and averages 
of the & equal to zero. The contribution from the collision part would vanish as one takes 
the limit because in the first collision term < Q(z — e) >= 1 and in the second collision term 
z is replaced by e due to the presence of the delta function. Upon simplification the velocity 
of the centre of mass would come out as 

v ™ = -y^T) < z > < 13) 

which is exactly the same as the one we had got in the over-damped calculations. Note that 
finding < V z > is simple in the present case because we know the velocity distributions of 
v\ and t>2 and < V z >=< v\ — v 2 > can be checked to be vanishing. Note that < V z > 
cannot be found out from the Eq.10 because it does not take into account the part of V z 
that comes during the collision happens and that is why < Q(z — e)V z > has a non- vanishing 
value whereas < V z > be vanishing in equilibrium. Its not difficult to see that, had we taken 
the collision contribution to the V z considering the presence of damping during collision, the 
collision terms contributing to this part of V z would have a factor 1/(1 — Q(z — e)) instead 



of l/Q(z — e) as is there for V cm . Now, the collision part would become non-vanishing unlike 
the case of its contribution to the V cm . Before we go for a generalization of our model to 
the soft collision case, let us first compare our result with that obtained from the direct 
numerical simulation of the full model comprising of Eq.1-3. 

Fig.l shows four plots of the ln(|V^ m |) against ln(a) at four different values of k^T . 
All these plots are straight lines with slope almost equal to 0.5. The match between the 

n 

theoretical and numerical values are quite good. In the previous communication [1[ where 
we had mentioned that there is some order of magnitude mismatch between the theoretical 
and numerical value of V cm that was because of not taking proper scaling of the noise term 
with time interval. In the over-damped case, Ax = yj2ksT /Yy/ At where T is the damping) 
is basically the statement of Fick's Law of diffusion. So, any numerical value put in the 
simulation for the ksT can be directly taken for analytic calculations. But, for the full 
(containing inertial term) model simulation we use Am = ^J2YkBT^fAi where u is velocity. 
We have to convert this expression to the equivalent form of the Fick's Law and then calculate 
the kftT and that would make k%T = At 3 n where we have put the number n in the place 
of ksT in the numerical analysis. In this rescaling At 2 comes from the conversion of Am to 
Ax and the rest comes because of the gamma appearing in the numerator of the noise term 
mentioned above. Note that in the numerical simulation to implement hardcore collision we 
simply exchange the position and velocity of the particles following any overshooting and 
that we can do because of considering the same mass for the particles. The time step in our 
simulation has been taken At = 0.01. 

Having understood the model with hardcore collisions, let us generalize it for soft colli- 
sions. At the core of the symmetry breaking in the model is the consideration of two different 
damping, one during collision and the other when the particles are not in touch with each 
other but are interacting harmonically. Let us consider that the particles have a damping 
T' when they collide which is different from that they have on flight (namely (1 — /3) and 
1). This may happen because when particles are actually in touch with each other their 
surroundings are different from that when they are on flight. This T' can also be different 
for different particles of the dimer, but, for the sake of simplicity one may consider them 
to be equal and that would not change the qualitative result. Now considering the collision 
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FIG. 1. Plot of ln(\v cm \) vs ln(a) at /3 = 0.9 for values of k B T = 5 x 1(T 7 , k B T = 5 x 1(T 9 , 
k B T = 5 x 10 -11 , and k B T = 5 x 10~ 13 respectively from top to bottom. The straight lines are 
plots of theoretical values and symbols correspond to numerical results. 



potential is soft and of width e one can rewrite the model as 



v x = B(Z - e)[ - (1 - /3) Vl + y/2(l - P)k B T Vl (t)] 

9 \e {z -e)^ + (l-e(z-e))(a-^) 



dz 



[I - Q(z - e))[T'v 1 - ^2V'k B Tr ll {t)\ 



v 2 = Q(Z - e)[ - v 2 + y/2k B Trj 2 (t)} 



d_ 

dz 



Q(z-e)^- + (l-e(z-e)){a-^) 
_(1 _ Q (z - e))[T'v 2 - y/2rk B T V2 (t)] 



Z = Vi - v 2 . 



(14) 



(15) 



Notice that, the temperature remains the same (T) irrespective of collision or harmonic 
interaction is taking place. This model can be analyzed similarly as the one shown above 
and would result in a centre of mass motion of the system over a range of values of V 
and P even showing conditions for the directional reversal of V cm . We will show a detailed 



analytic and numerical treatment of the above model in future. Imagine the fact that, the 
dimer and the bath composite system is thermally insulated from the surroundings. Now, 
also imagine that at collision some other internal degrees of freedom of the particles are 
excited (particles are not dimensionless particles) and that results in an inelastic collision of 
the particles. Eventually, these internal coordinates would relax leaving the energy back to 
the bath and consider that this process happens pretty quickly. Now, this isolated system 
must reach equilibrium over time since its an isolated system from the surroundings, but, 
intermittent stealing of energy by the internal degrees of freedoms of the particles would 
result in a kinetic energy oscillation of the dimer. This kinetic energy oscillation can be 
imagined as an oscillation of the effective temperature of the particles at shorter time scales 
smaller than the one over which the average temperature for the equilibrium is defined. One 
can even tune the time period of this effective temperature oscillation of the dimer by tuning 
the force constant a because it sets the time period of oscillation of the dimer on average. 
Not only that, this time period can also be large. This is a typical so-called non-equilibrium 
state (on average steady oscillation) which would be coexisting with equilibrium because 
the isolated system must reach equilibrium at large time. So, this would be a possible way 
of getting a non-equilibrium order within the time scale over which equilibrium is defined 
and its not all uncorrelated fluctuations within the equilibrium time scale. There will be 
associated positive and negative productions of entropy (rather entropy exchange between 
bath and the particles) such that the total entropy production rate over a cycle remains zero 
to maintain equilibrium at larger time scales. So, this system is an entropy oscillator which 
is consistent with equilibrium defined over a somewhat larger scale. 
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